Block boundary artifact reduction for block-based image compression

ABSTRACT

Block boundary artifact reduction in decompressed digital images is accomplished by filtering the intensity of pixels in the vicinity of the block boundary. The filter utilizes filter coefficients selected from tables on the basis of the distribution of a scalar quantity describing the pixels neighboring the pixel to which the intensity adjustment is to be applied. The intensities of pixels at the boundary and one pixel removed from the boundary are adjusted. If interpolation pixels in each of neighboring blocks are of relatively constant intensity, the intensities of additional pixels, more remote from the boundary, are adjusted. Filter coefficients can be selected from different arrays for pixels on the boundary or removed from the boundary or if the pixel intensity adjustment is based on the intensities of pixels in a horizontal row or vertical column of interpolation pixels.

BACKGROUND OF THE INVENTION

The present invention relates to digital image reproduction and more particularly to a technique for reducing block boundary artifacts in digital images which have been subjected to block based image compression and decompression.

In various applications, including high definition television and video conferencing systems, it is desirable to store and transmit images in a digital form. The amount of data generated by converting images from analog to digital form is so great that digital transmission and storage of images would be impractical if the data could not be compressed to require less bandwidth and data storage capacity. As a result, a number of techniques have been developed to compress the digital data representing an image and a video sequence.

Several image and video compression techniques process digital image data by first dividing the image into 2-dimensional blocks. A common block size for block-based compression is an area eight pixels by eight pixels in size. Images and video may be processed in 8×8 blocks by algorithms specified by the ISO “JPEG” specification for still images, the CCITT “H.261” recommendation for video conferencing, and the ISO “MPEG” specification for video compression. After the image is divided into blocks, compression of the data proceeds on a block-by-block basis. There are several methods for processing the intrablock data and more than one method may be used. A common compression technique includes a transformation step where a transform, such as a discrete cosine transform (DCT), is applied to the intrablock image data to generate transform coefficients for each block of picture elements or pixels. The transformation is often followed by a quantization step where the transformed intrablock image data is replaced by quantized data which usually requires fewer data bits. The quantized data may be processed further and then stored or transmitted as compressed image data. When the compressed image data is received, the receiver reverses the process, decompressing the data and reconstructing the image.

The image compression methods with the highest compression efficiency are so called “lossy” methods. An image often contains considerable information, such as small detail, that contributes very little to a real observer's appreciation of the reconstructed image. When transformed this information appears as intermediate or higher frequency components which can be discarded following transformation. In addition, some degradation of the image may be acceptable in certain circumstances as a trade-off to reduce the amount of data which must be transmitted or stored. Discarding or losing “excess” information during the quantization of the transformed data is the core of the data compression effectiveness of lossy compression methods. On the other hand, much of the information discarded in block-by-block compression is data relevant to interblock relationships. Processing image data on a block-by-block basis and discarding information related to interblock relationships often results in a noticeable change in pixel intensity across the boundary between blocks of the reconstructed image. The “blocky” appearance of a decompressed image has been consistently recognized as a problem when using lossy compression methods for digital image processing.

Several strategies have been developed to reduce the block boundary artifacts in a decompressed digital image. For example, the H.261 compression—decompression technique incorporates an optional low pass filter that can be applied to all pixels of the image except those at the edge of the image. Applying the filter to all pixels of the image is time consuming and reduces the apparent image resolution. Another technique is the projection on convex sets (POCS) method which attempts to estimate the characteristics of the pre-compression image from the compressed data together with a model of the statistical characteristics of the image before compression. This method is successful in producing high quality decompressed images and eliminating block boundary artifacts. However, the computational complexity of the method makes it impractical for real time video applications.

Another strategy to reduce block boundary artifacts involves retaining some of the interblock correlation data when the image is compressed. Lossy compression commonly involves a transformation step where spatial domain image data is converted to a frequency domain representation by application of a mathematical transform, such as a discrete cosine transform (DCT). The DCT coefficients have been used to determine whether to apply a filter to the image data and, if so, to determine which filter to apply. The interblock correlation data can either be embedded in the transform coefficients or transmitted separately. However, the DCT coefficients which are not incorporated in the compressed data must be transmitted and stored separately from the image data reducing the efficiency of the data compression. Further, integration of this additional data in existing standards, such as the JPEG standard, is difficult or impossible.

Komuro et al., U.S. Pat. No. 5,675,666, disclose a technique for preprocessing image data before compression to compensate for the “block effect” artifact that would normally appear upon decompression of the image data. Error correction data is developed by a test compression and decompression of the image followed by computation of the error correction data from the differences between the input image data and the data for an output test image. The error correction data is added to the input image data so that when the image is compressed again and decompressed, subsequently, adjustments to reduce block boundary artifacts will have been made in the image data for pixels at the block boundaries. Preprocessing can be effective in reducing the boundary artifacts but requires access to the image content before compression.

Images may originate from many sources which do not have access to preprocessing capability. As a result, techniques of postprocessing the image to reduce block boundary artifacts after decompression in order to restore or enhance image quality have been proposed. One postprocessing technique of boundary artifact reduction is disclosed by Kim, U.S. Pat. No. 5,694,492. In this technique the compressed image data is received, decompressed and stored. A target pixel is identified and its data is filtered repeatedly until the absolute difference between the filtered and unfiltered data exceeds a predefined threshold. The process is repeated for each pixel in the image. This technique relies on repetitive filtering of data with a low pass filter having fixed coefficients. The iterative nature of the process is time consuming and computationally expensive. In addition, the fixed filter coefficients do not respond to local signal variations in the image resulting in less than optimal treatment of some discontinuities.

Chan et al., in a paper entitled A PRACTICAL POSTPROCESSING TECHNIQUE FOR REAL-TIME BLOCK-BASED CODING SYSTEM, IEEE Transactions on Circuits and Systems for Video Technology, Vol. 8, No. 1 (February 1998) describe a postprocessing boundary artifact reduction technique based on a data-dependent, spatially-variant filter which operates on the pixels on both sides of the block boundaries. Filter coefficients are based on a vector quantization of pixel intensity values of several neighboring pixels on each side of the block boundary in the vicinity of the boundary pixel to which intensity adjustment will be applied. While this technique adjusts the filter in response to the local signal conditions, vector quantization is slow and computationally intensive. In essence, vector quantization is the comparison of a vector of pixel intensity values to a set of n dimensional vectors to determine the most representative vector. Further, the technique adjusts the intensity of only those pixels adjoining the boundary. This limits the available smoothing of the intensity change to reduce the boundary artifact. The limited depth of the adjustment zone is particularly apparent at boundaries between blocks of relatively constant intrablock intensity. However, adjusting more than a single row of pixels on each side of the boundary would substantially increase the time and computational resources necessary to perform artifact reduction with this vector-based method.

What is desired, therefore, is a post processing technique that requires no transmission or storage of transform coefficients, can be performed rapidly enough to be useful with digital video, is computationally conservative, and adjusts for local signal variations within the image to effectively reduce block boundary artifacts under a variety of circumstances.

SUMMARY OF THE INVENTION

The present invention overcomes the aforementioned drawbacks of the prior art by providing a method for reducing an artifact at a block boundary in a decompressed digital image including the steps of identifying a plurality of contiguous interpolation pixels including at least one pixel on each side of the block boundary; identifying one interpolation pixel on one side of the block boundary as a subject pixel; determining the distribution of a scalar characteristic of the interpolation pixels; selecting a set of filter coefficients correlated to the distribution of the scalar characteristic; modifying the intensity of the subject pixel by applying a filter producing the modified intensity as a function of the selected filter coefficients. In a preferred embodiment, the intensity of both a border pixel adjoining the block boundary and a neighboring pixel, adjoining the border pixel but one pixel removed from the boundary, are adjusted. Adjusting the intensity of pixels in a zone two pixels deep reduces the abruptness of the transition between blocks and improves block boundary artifact reduction by the linear filter. Basing the filter on a scalar characteristic of the pixel, such as its intensity, produces a method which is computationally conservative and can be performed rapidly enough to be used with video.

In a first aspect of the invention, different filter coefficients are applied if the subject pixel is a border pixel adjoining the block boundary or a neighboring pixel, one pixel removed from the block boundary. Providing a less aggressive adjustment to pixels remote from the block border further reduces the abruptness of the transition. In a second aspect of the invention, filter coefficients are selected from different tables of coefficients if the intensity adjustment is based on interpolation pixels in a horizontal row or a vertical column. Images characteristically have a different texture in the horizontal and vertical directions. In a third aspect of the invention, the intensities of additional pixels, more remote from the block boundary, are adjusted if the intensities of the interpolation pixels are substantially constant in the neighboring blocks. Increasing the depth of the adjustment zone along the block boundary is important to reducing the block boundary artifact between two relatively large areas of constant intensity.

The foregoing and other objectives, features and advantages of the invention will be more readily understood upon consideration of the following detailed description of the invention, taken in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic illustration of an exemplary video processing block.

FIG. 2 is a flow chart illustrating a technique for selecting filter coefficients for a pixel intensity filter.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

Many digital image compression schemes, such as the ISO “JPEG” specification for still images and the ISO “MPEG” specification for video, rely, at least in part, on block-based image compression. Referring to FIG. 1, the image is divided into a series of two dimensional areas or blocks 2, indicated by a bracket, and processing proceeds on a block-by-block basis. A common size processing block is eight pixels 4 in height by eight pixels 4 in width. Following transmission or storage, the image is reconstructed by reversing the process to decompress the image on a block-by-block basis. Data compression efficiency can be enhanced by discarding some information related to the image which is considered unnecessary for a viewer's appreciation of the reconstructed image. The discarded information is generally related to intermediate or higher frequency components of the data which describe minor details of the image. With block-by-block compression techniques, a part of the discarded information describes the relationship between a block 2 and its neighbors. As a result, the intensity of the pixels in the decompressed image can change dramatically across the block boundaries 6. These artifacts of the block-by-block compression can result in visible blocks in the decompressed image. In the preferred embodiment, a spatially-variant linear filter is used to modify the intensities of pixels 4 near the block boundary 6 to reduce block boundary artifacts and smooth the transition between blocks. The coefficients of the filter are selected to correspond to the distribution of intensities of a group of pixels contiguous to the pixel to which the modification is to be applied. As a result, the filtering applied will be appropriate to the type and magnitude of the block boundary artifact present at that particular point in the image.

The linear filter of the preferred embodiment is of the general form: ĥ _(j) =c _(ai) h _(i) +c _(aj) h _(j) +c _(ak) h _(k) +c _(al) h _(l) where:

-   -   i, j, k, and l are a plurality of contiguous interpolation         pixels crossing a block boundary in a horizontal row or vertical         column;

-   ĥ_(j)=the adjusted intensity of the subject pixel, j;

-   h_(j)=the unmodified intensity of the subject pixel, j, in the     decompressed image; and

-   c_(ax)=a filter coefficient applicable to the intensity of the     specified pixel.     The spatially variant linear filter is applied to adjust the     intensity of pixels in the block border rows 8 on both sides of and     adjoining to the block boundary 6. Improved artifact reduction is     achieved by also applying a spatially variant linear filter to     adjust the intensity of pixels in the neighboring rows 10, one pixel     removed from the block boundary 6. Increasing the depth of the     intensity adjustment zone further reduces the abruptness of the     intensity transition and further smooths the transition between     blocks.

It is to be understood that the particular filter selected may be linear or non-linear. The filter may be multi-dimensional although a one dimensional filter is preferred. Also, the filter may be an arbitrary width so long as it preferably includes the respective pixel being examined.

To apply the technique of block boundary artifact reduction of the preferred embodiment a group of contiguous interpolation pixels in a horizontal row 12, 14, 16, 18, 20, 22, 24, and 26 or a vertical column 28 of pixels (having intensities v₀ through v₇) is identified in a decompressed image, see FIG. 1. The interpolation pixels could be non-contiguous and not in a horizontal row or vertical column, if desired. One of the interpolation pixels in either the block boundary row 8 of pixels or the neighboring row 10 of pixels is identified as the subject pixel to which intensity adjustment will be applied. For example, the pixel 16 having intensity h₂ located in the block boundary row 8 is subject to intensity adjustment to reduce block boundary artifacts. The intensities of the subject pixel 16 and three contiguous neighboring pixels of intensities h₁ 14, h₃ 18, and h₄ 20 are utilized in adjusting intensity of the subject pixel 16. In general, the intensities of the subject pixel 16, its contiguous neighbor 14 in the direction away from the block boundary 6 and its two contiguous neighboring pixels 18 and 20 in the direction of the block boundary identified 40 and incorporated into the intensity adjustment filter, see FIG. 2. The filter appropriate for adjusting the intensity of the subject pixel 16 with intensity h₂ would include the intensities of pixels 14, 16, 18 and 20 and be of the form: ĥ ₂ =c _(a1) h ₁ +c _(a2) h ₂ +c _(a3) h ₃ +c _(a4) h ₄

Likewise the filters appropriate for adjusting the intensities of the pixels 18, 20, and 22 would be, respectively: ĥ ₃ =c _(b1) h ₂ +c _(b2) h ₃ +c _(b3) h ₄ +c _(b4) h ₅ ĥ ₄ =c _(c1) h ₅ +c _(c2) h ₄ +c _(c3) h ₃ +c _(c4) h ₂ ĥ ₅ =c _(d1) h ₆ +c _(d2) h ₅ +c _(d3) h ₄ +c _(d4) h ₃

The intensities of the subject pixel 16 and its three contiguous neighbors 14, 18, and 20 are also used to select the appropriate filter coefficients, c_(ai), from a table of filter coefficients. In the case of the example subject pixel 16, the mean of the pixel intensities h₁ 14, h₂ 16, h₃ 18, and h₄ 20 is calculated 42. An intensity differential, the difference between the intensity of each of the four interpolation pixels and the mean intensity, is then calculated 44. By definition, if the intensities are not equal, the mean of the intensities must be greater than some of the individual intensities and less than others. Therefore, the intensity differentials will be both positive and negative.

It is to be understood that any other suitable statistical measure may be used to determine intensity variations.

The calculated intensity differentials are quantized 46 by sorting into four bins; negative and less than the negative of a threshold (−t) 48, negative but greater than or equal to the negative of the threshold; zero or positive but less than or equal to the positive of the threshold (t) 50; and positive and greater than the positive threshold. Experimentation with sample images indicates that a threshold intensity differential (t) between four and eight is appropriate with six being optimal. Each bin is represented by a two-bit number 51. The four, 2-bit, quantized intensity differentials are concatenated 52 to form a single 8-bit binary index value 54. The 8-bit index value 54 represents the intensity distribution of the four interpolation pixels for that subject pixel. The 8-bit index value 54 is used to select 56 appropriate filter coefficients 57 from a table 58 of index values 54 and corresponding filter coefficients 57. Since the index value 54 represents the intensity distribution of pixels in the vicinity of the subject pixel to which adjustment will be applied, the filter coefficients 57 selected for use in modifying the intensity of the subject pixel will be an optimal set of coefficients for the block boundary artifact at that point in the image. Since intensity differentials are both positive and negative, 81 combinations of index values can be eliminated as “not possible,” which reduces the necessary size of the table.

It is to be understood that the number of “bins” into which the intensity differentials are mapped may be varied as desired. In addition, the values of the boundaries between the different “bins” may be selected as desired. Further, bin and bin boundaries need not be symmetrical or symmetrical about zero. Also, the technique of concatenation to obtain an index value may be replaced by any suitable technique to define the result of quantizing the intensity differentials. One technique for establishing appropriate values for filter coefficients, is to compress and decompress a series of test images with a lossy compression-decompression technique or codec to be used in the application. After the images have been compressed and then decompressed, the block boundaries 6 of the reconstructed image can be analyzed by comparing data related to the original image with that for the points at which artifact reduction might be applied in the “blocky” decompressed image. A technique such as the least-squares technique can be used to compute an optimal set of filter coefficients 57 for each set of data corresponding to a quantized index value derived from the intensity distribution of interpolation pixels crossing the block boundary. This can be accomplished by inverting the sense of the filter equations, above, using the intensities ĥ_(i) and h_(i) from actual pixel data and solving for the appropriate filter coefficients, c_(ai), c_(bi), c_(ci), and c_(di).

Another technique of establishing the values of filter coefficients 57 would be to calculate the coefficients. Each index value implies an intensity distribution of the corresponding interpolation pixels. The intensity distribution implied by the index value can be used to imply intensities of individual interpolation pixels for the purposes of calculating the appropriate filter coefficients. Any other suitable technique may likewise be used.

As illustrated in FIG. 1, border pixels 8 and neighboring pixels 10 may be part of a horizontal row 30 or a vertical column 28 of interpolation pixels. In addition, there are a number of pixels located in block border rows 8 or neighboring pixel rows 10 for both a horizontal and a vertical block border 6. The intensity of these pixels are filtered in both horizontal and vertical directions. The present inventor came to the realization that the intensity smoothing at the block boundary 6 may be improved by increasing the depth of the adjustment zone, that is, by adjusting the intensity of both the border pixel 8 and the pixels in the row one pixel removed from the border, the neighboring pixel 10. While the same filter coefficients can be used, the present inventor came to the realization that boundary artifact reduction can be further enhanced by selecting filter coefficients from different tables providing a more aggressive intensity adjustment for border pixels 8 than for neighboring pixels 10 which are further from the boundary. Likewise, filter coefficients may be selected from the same array of filter coefficients for pixels which will be adjusted on the basis of belonging to interpolation pixels in a horizontal row 30 or a vertical column 28. The present inventor's experimentation has shown that mixing the data obtained from horizontal and vertical sampling of test data to obtain a common set of coefficients affords excellent image quality and minimizes the data storage requirements. However, the human vision system anticipates different textures in images and nature in the horizontal and vertical directions and the filter coefficients used when processing pixels in a horizontal row 30 of interpolation pixels can be varied from those used when processing pixels in a vertical column 28 of interpolation pixels. To process pixels in the vicinity of the corners of the block, filtering can be performed in one direction, for example horizontal, followed by filtering in the other direction, for example vertical, using the filtered intensities from the first operation as input for the processing in the second direction.

When neighboring blocks of the decompressed image are of relatively constant intensity additional filtering at the common block boundary 6 is desirable to visually improve the block boundary transition. The present inventor came to the realization that increasing the depth of the adjustment zone by adjusting pixels further removed from the boundary improves the transition between the relatively large areas of uniform intensity. Referring to FIG. 1, if the group of interpolation pixels 12, 14, 16, and 18 and the group of interpolation pixels 20, 22, 24, and 26 in a neighboring block are each of substantially constant intensity adjusting the intensity of the pixels 12 and 14, three and four pixels removed from the boundary 6, in addition to the border pixel 18 and the neighboring pixel 16 enhances the reduction of the boundary artifact. Likewise, in addition to the border pixel 20 and neighboring pixel 22, pixels 24 and 26 would be adjusted in the neighboring block. If the intensity of the interpolation pixels are substantially equal in each block, the following filter equations are appropriate for pixels 12, 14, 16, and 18, respectively: ĥ ₀=(15h ₃ +h ₄)/16 ĥ ₁=(13h ₃+3h ₄)/16 ĥ ₂=(11h ₃+5h ₄)/16 ĥ ₃=(9h ₃+7h ₄)/16

The terms and expressions that have been employed in the foregoing specification are used as terms of description and not of limitation, and there is no intention, in the use of such terms and expressions, of excluding equivalents of the features shown and described or portions thereof, it being recognized that the scope of the invention is defined and limited only by the claims that follow. 

1. A method of reducing an artifact at a block boundary in a decompressed image comprising the steps of: (a) identifying a plurality of interpolation pixels including at least one pixel on each of a first side and a second side of said block boundary; (b) identifying one said interpolation pixel on said first side of said block boundary as a subject pixel; (c) determining a distribution of a scalar characteristic of said interpolation pixels, wherein said scalar characteristic is a non-quantized characteristic; (d) selecting a plurality of filter coefficients correlated to a quantization of said distribution of said scalar characteristic; (e) modifying an intensity of said subject pixel by applying a filter producing said modified intensity of said subject pixel as a function of said filter coefficients; and (f) selecting a second plurality of filter coefficients correlated to a distribution of said scalar characteristic of said first pixel, said subject pixel, a second pixel on said second side of said block boundary and adjoining said subject pixel, and a third pixel adjoining a side of said first pixel remote from said block boundary; and (g) modifying an intensity of a first pixel by applying a filter producing said modified intensity of said first pixel as a function of said second plurality of filter coefficients.
 2. A method of reducing an artifact at a block boundary in a decompressed digital image comprising the steps of: (a) identifying a plurality of interpolation pixels including at least one pixel on each of a first side and a second side of said block boundary; (b) identifying one said interpolation pixel adjoining said first side of said block boundary as a border pixel; (c) identifying one said interpolation pixel adjoining a side of said border pixel remote from said block boundary as a neighboring pixel; (d) modifying an intensity of said border pixel by application of a first filter producing said modified intensity as a function of filter coefficients correlated to a quantized distribution of non-quantized intensities of a first plurality of said interpolation pixels; and (e) modifying an intensity of said neighboring pixel by application of a second filter producing said modified intensity as a function of filter coefficients correlated to a quantized distribution of non-quantized intensities of a second plurality of said interpolation pixels.
 3. The method of claim 2 wherein said first plurality of interpolation pixels comprises: (a) said border pixel; (b) said neighboring pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) a fourth interpolation pixel adjoining a side of said third interpolation pixel remote from said block boundary.
 4. The method of claim 2 wherein said second plurality of interpolation pixels comprises: (a) said neighboring pixel; (b) said border pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) a fourth interpolation pixel adjoining a side of said neighboring pixel remote from said block boundary.
 5. The method of claim 2 wherein said interpolation pixels of said second plurality are said interpolation pixels of said first plurality.
 6. The method of claim 2 wherein said second plurality of interpolation pixels includes at least one said interpolation pixel not included in said first plurality of interpolation pixels.
 7. The method of claim 2 wherein said second filter produces said modified intensity of said neighboring pixel as a function of filter coefficients correlated to a distribution of intensities of said first plurality of said interpolation pixels.
 8. A method of reducing an artifact at a block boundary in a decompressed digital image comprising the steps of: (a) identifying a plurality of interpolation pixels including at least one pixel on each of a first side and a second side of said block boundary; (b) identifying one said interpolation pixel adjoining said first side of said block boundary as a border pixel; (c) identifying one said interpolation pixel adjoining a side of said border pixel remote from said block boundary as a neighboring pixel; (d) selecting a plurality of border pixel filter coefficients correlated to a non-quantized distribution of intensities of a first plurality of said interpolation pixels from a first array of filter coefficients; (e) selecting a plurality of neighboring pixel filter coefficients correlated to a non-quantized distribution of intensities of a second plurality of said interpolation pixels from a second array of filter coefficients; (f) modifying said intensity of said border pixel by application of a first filter producing said modified intensity as a function of said border pixel filter coefficients; and (g) modifying said intensity of said neighboring pixel by application of a second filter producing said modified intensity as a function of said neighboring pixel filter coefficients.
 9. The method of claim 8 wherein said first plurality of interpolation pixels comprises: (a) said border pixel; b) said neighboring pixel; c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and d) a fourth interpolation pixel adjoining a side of said third interpolation pixel remote from said block boundary.
 10. The method of claim 8 wherein said second plurality of interpolation pixels comprises: (a) said neighboring pixel; (b) said border pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) a fourth interpolation pixel adjoining a side of said neighboring pixel remote from said block boundary.
 11. The method of claim 8 wherein said interpolation pixels of said second plurality are said interpolation pixels of said first plurality.
 12. The method of claim 8 wherein said second plurality of interpolation pixels includes at least one said interpolation pixel not included in said first plurality of interpolation pixels.
 13. The method of claim 8 wherein said filter coefficients of said second array are said filter coefficients of said first array.
 14. The method of claim 8 wherein said second plurality of filter coefficients includes at least one said filter coefficient not included in said first plurality of filter coefficients.
 15. The method of claim 8 wherein said neighboring pixel filter coefficients are equal to said border pixel filter coefficients.
 16. The method of claim 8 wherein said neighboring pixel filter coefficients include at least one filter coefficient which is not equal to at least one of said border pixel filter coefficients.
 17. A method of reducing an artifact at a block boundary in a decompressed digital image comprising the steps of: (a) identifying a plurality of interpolation pixels including a first plurality of at least three said interpolation pixels on a first side of said block boundary and a second plurality of at least three said interpolation pixels on a second side of said block boundary; (b) identifying one said interpolation pixel adjoining said first side of said block boundary as a border pixel; (c) identifying one said interpolation pixel adjoining a side of said border pixel remote from said block boundary as a neighboring pixel; (d) identifying a said interpolation pixel adjoining a side of said neighboring pixel remote from said block boundary as a remote neighboring pixel; (e) comparing a first plurality of non-quantized intensities of said first plurality of interpolation pixels to each other; (f) comparing a second plurality of non-quantized intensities of said second plurality of interpolation pixels to each other; (g) modifying said intensities of said border pixel and said neighboring pixel if at least one of said first plurality of intensities are not substantially equal and said second plurality intensities are not substantially equal; and (h) modifying said intensities of said border pixel, said neighboring pixel, and said remote neighboring pixel if said first plurality of intensities are substantially equal and said second plurality of intensities are substantially equal.
 18. The method of claim 17 wherein said intensities of said border pixel and said neighboring pixel are modified if said intensities of at least one of said first plurality of intensities or said second plurality intensities are not substantially equal by the steps of: (a) selecting a plurality of border pixel filter coefficients correlated to a distribution of intensities of a third plurality of interpolation pixels; (b) modifying said intensity of said border pixel by application of a first filter producing said modified intensity as a function of said border pixel filter coefficients; (c) selecting a set of neighboring pixel filter coefficients correlated to a distribution of intensities of a fourth plurality of interpolation pixels; and (d) modifying said intensity of said neighboring pixel by application of a second filter producing said modified intensity as a function of said neighboring pixel filter coefficients.
 19. The method of claim 18 wherein said neighboring pixel filter coefficients are equal to said border pixel filter coefficients.
 20. The method of claim 18 wherein said neighboring pixel filter coefficients include at least one filter coefficient which is not equal to at least one of said border pixel filter coefficients.
 21. The method of claim 18 wherein said third plurality of interpolation pixels comprises: (a) said border pixel; (b) said neighboring pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) a fourth interpolation pixel adjoining a side of said third interpolation pixel remote from said block boundary.
 22. The method of claim 18 wherein said fourth plurality of interpolation pixels comprises: (a) said neighboring pixel; (b) said border pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) said remote neighboring pixel.
 23. The method of claim 17 wherein if said first plurality of intensities are substantially equal and said second plurality of intensities are substantially equal said intensities of said border pixel, said neighboring pixel, and said remote neighboring pixel are modified by the steps of: (a) modifying said intensity of said border pixel by application of a first filter producing said modified intensity as a function of said intensity of said border pixel and an intensity of one said interpolation pixel on said second side of said block boundary and adjoining said border pixel. (b) modifying said intensity of said neighboring pixel by application of a second filter producing said modified intensity as a function of said intensity of said border pixel and an intensity of one said interpolation pixel on said second side of said block boundary and adjoining said border pixel. (c) modifying said intensity of said remote neighboring pixel by application of a third filter producing said modified intensity as a function of said intensity of said border pixel and an intensity of one said interpolation pixel on said second side of said block boundary and adjoining said border pixel.
 24. A method of reducing an artifact at a block boundary in a decompressed digital image comprising the steps of: (a) identifying a plurality of interpolation pixels including at least one pixel on each of a first side and a second side of said block boundary (b) identifying one said interpolation pixel adjoining said first side of said block boundary as a border pixel; (c) identifying one said interpolation pixel adjoining a side of said border pixel remote from said block boundary as a neighboring pixel; (d) selecting a plurality of border pixel filter coefficients correlated to a distribution of intensities of a first plurality of said interpolation pixels from a first array of filter coefficients if said interpolation pixels are arrayed horizontally in said image; (e) selecting a plurality of neighboring pixel filter coefficients correlated to a distribution of intensities of a second plurality of said interpolation pixels from a second array of filter coefficients if said interpolation pixels are arrayed horizontally in said image; (f) selecting said plurality of border pixel filter coefficients from a third array of filter coefficients if said plurality of interpolation pixels are not arrayed horizontally in said image; (g) selecting said plurality of neighboring pixel filter coefficients from a fourth array of filter coefficients if said plurality of interpolation pixels are not arrayed horizontally in said image; (h) modifying an intensity of said border pixel by application of a first filter producing said modified intensity as a function of said border pixel filter coefficients; and (i) modifying an intensity of said neighboring pixel by application of a second filter producing said modified intensity as a function of said neighboring pixel filter coefficients.
 25. The method of claim 24 wherein said filter coefficients of said second array are equal to said filter coefficients of said first array and said filter coefficients of said third array are equal to said filter coefficients of said fourth array.
 26. The method of claim 24 wherein said filter coefficients of said third array are equal to said filter coefficients of said first array and said filter coefficients of said fourth array are equal to said filter coefficients of said second array.
 27. The method of claim 24 wherein said filter coefficients of said second array include at least one filter coefficient which is not equal to at least one of said filter coefficients of said first array and said filter coefficients of said third array include at least one filter coefficient which is not equal to at least one of said filter coefficients of said fourth array.
 28. The method of claim 24 wherein said filter coefficients of said third array include at least one filter coefficient which is not equal to at least one of said filter coefficients of said first array and said filter coefficients of said second array include at least one is-filter coefficient which is not equal to at least one of said filter coefficients of said fourth array.
 29. The method of claim 24 wherein said first plurality of interpolation pixels comprises: (a) said border pixel; (b) said neighboring pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and (d) a fourth interpolation pixel adjoining a side of said third interpolation pixel remote from said block boundary.
 30. The method of claim 24 wherein said second plurality of interpolation pixels comprises: (a) said neighboring pixel; (b) said border pixel; (c) a third interpolation pixel adjoining a side of said border pixel nearest said block boundary; and d) a fourth interpolation pixel adjoining a side of said neighboring pixel remote from said block boundary. 